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Abstract 

Consider a multiple component system connected in parallel. In this system, 
as components fail one by one, the total load or traffic applied to the system is 
redistributed among the remaining surviving components, which is commonly 
referred to a s load-sharing. 



Recently 



Singh et al. 



(J2008I ) proposed different 



Kim and Kvaml (J2004I ) and 
load-sharing models and developed parametric inference for the these models. 
However, their parametric estimates are calculated using iterative numerical 
methods. In this note, we provide the general closed-form MLEs for the two 
load-sharing models provided by them. 

Keywords: Reliability, load-sharing, maximum likelihood estimate 
(MLE), closed- form solution. 



1 Introduction 

Most research work involving load sharing models has mainly focused on the char- 
acterization of system reliability under a known load-sharing rule and parameters. 
The parameter estimation of the load-sharing rule has not yet been fully developed. 
Recently, parametr ic inference for r eliabil ity u nder the equal loa d-sharing rule has 



been considered by 



Kim and Kvaml (120041 ) and 



Singh et al. 



(120081 ). They solved the 



likelihood estimating equations to find the maximum likelihood estimators (MLEs) 
of the load-sharing parameters. 

However, they provide no general closed form solutions for the MLEs, but instead 
use iterative numerical methods to calculate their estimates. It is well known that 
there are some problematic issues associated with iterative numerical methods such 
as stability and convergence. 



In this note, we provide the general c 



models provided by 



osed- 
Kim and Kvaml (J2004J) and 



orm MLEs f or th e two load-sharing 



Singh et al, 



tosh . 



2 Kim-Kvam load-sharing model 



Consi der /c-component system connected in parallel. Following iKim and Kvam 



( 120041 ). we assume the following: 



(i) A system is made up of k components whose lifetimes are independent and have 



identical exponential distributions with initial failure rate 9. 

(ii) After the first component fails, the failure rates of the remaining k — 1 com- 
ponents change to Ai$ where Ai > 0. After the next component failure, the 
failure rates of the surviving k — 2 components change to ^9. Then, after the 
next component failure, the the failure rates of the surviving k — 3 components 
change to \%9 and so on and so forth. 

(iii) There are n repeated measurements of independent systems. That is, we have 
a random sample of independent systems of size n. 



Let X im denote the lifetime of the m-th component in the 2-th parallel system where 
i = 1,2, ...,n and m = 1,2, ... ,k. For notational convenience, we can re-index 
the lifetimes such that Xn < Xj2 < ■ • • < X^. Then the time spacing between the 
(j — l)-th failure and j-th failure for the i-th system is T^ = Xjj — Xjj_i with X i0 = 0. 



As shown in 



Kim and Kvaml (J2004f ). the likelihood function for the i-th system is 



L i (9,A) = (k\)9 k -[f[\^ 1 

3=1 



■ exp 






i.i 



where A = 1 and A = (Ai, . . . , A^_i). It is immediate that the likelihood function for 
a random sample of size n is given by 

k n k 



L(9,A) = (k\) n 9 nk - [ JA^_! 



• exp 






ij 



Taking the logarithm of ([TJ, differentiating with respect to 9, Ai, . . . , Afe_i, denoting 
partial derivative of log L with respect to 9 as Iq = d log L/89 and partial derivative 



of logL with respect to Xj-i as £j_i = d\ogL/d\j-i, we obtain the log-likelihood 
estimating equations shown below: 

n k 

e 



i=i j=i 



and 



ej-i = ^-Q X> - j + 1)^- = o (3) 

V 1 7=? 

for j = 2, 3, . . . , k. For convenience, we denote t»j = J27=i *y- Then, (J2J) and (J3J) can 



be rewritten as 



^ = ! ^-E^-^ + 1 )Vi^ = (4) 



and 



71 

i-i = T 0(fc-j + 1)^ = 0. (5) 



It is immediate from solving flHj) for A,_i that we have 



n 
Xj - 1= 9(k-j + l)t./ 3 = 2,.. .,k. ( 6 ) 



Since A = 1, we rewrite (TjJ as 



y - fci.i - ^(k ~ J + yXj-iLj = 0. (7) 

i=2 



Substituting (jHj) into ((Zj) gives 



nfc , n(A; — 1) 



Solving (JS} for 9, we obtain the MLE of 9, denoted by 9, 



a n n 

kt»i , v— \ 



HI 



The MLEs of \j-\, denoted by Aj_i, are also easily obtained by substituting (jHJ) into 

n 



A '-'° (i -j + l) f~ = ^~" J = 2,3,-,*- 






3 Singh-Shar ma- Kumar load-sharing model 



Singh et al. 



(|2008|) 



Consider /c-component system connected in parallel. Following 
we assume the following: 



(i) A system is made up of k components whose lifetimes are independent and have 
exponential distributions with initial failure rate 9. 

[ii) After the first component fails, the failure rates of the remaining k — 1 com- 
ponents change to Ai$ where Ai > 0. After the next component failure, the 
failure rates of the surviving k — 2 components change to A26 1 . Then, after the 
next component failure, the the failure rates of the surviving k — 3 components 
change to \%9 and so on and so forth. 

After the failure of a certain number of components, say, after the s-th compo- 
nent failure (s > 2), the failure rates of the k — s remaining components change 



to \ s t6 (linearly increasing failure rate). In a similar manner, after the (s + 1)- 
th component failure, the failure rates of the k — s — 1 remaining components 
change to X s+ it8, and so on. Finally, after the last failure, the failure rate of 
the last component becomes \k-it9. 

(hi) There are n repeated measurements of independent systems. That is, we have 
a random sample of independent systems of size n. 



Again, let Xi m denote the lifetime of the m-th component in the 2-th parallel 
system where i = l,2,...,n and m = l,2,...,k. For notational convenience, we 
can re-index the lifetimes such that X a < X i2 < • • ■ < X^. Then the time spacing 
between the (j — l)-th failure and j-th failure for the i-th system is Tjj = Xij — Xjj_i 
with X iQ = 0. 



As is given in 



Singh et al. 



(120081 ). the likelihood function for the 2-th system is 



L t (M)=(^'[nvi]-[ n ^ 

3=1 j=s+l 



x exp 



o{ j> - ^ + l ) X 3-iUo + 7; E ( k " i + x ) A i-i* 



3=1 j=s+l 

where Ao = 1 and A = (Ai, . . . , A ft _i). It is immediate that the likelihood function for 
a random sample of size n is given by 

ft n ft 



ft n ft 

L(e,A)=(wre Bfc .[nAy_ 1 ].[n n ** 



3=1 1=1 j = S + ] 

n s 



x exp 



E { I> - •?' + *) Vi*« + oE( fc ^' + ^Vi* 



i=l 3=1 



3=s+l 



(10) 



Taking the logarithm of ffTOj) . differentiating with respect to 9, Ai, . . . , A^-i, de- 
noting partial derivative of logL with respect to 9 as £ g = d log L/89 and partial 
derivative of logL with respect to Aj_i as £^\ = <91ogL/<9Aj_i, we obtain the log- 
likelihood estimating equations shown below: 

, n s _ k 

l » = t -E {!> -i+ ^Vi^ + 2^ (fc_J + i)V4) = o» (ii) 

i=l j=l J=s+1 

n 

ij-i=--^- -9j2(k-j + l)t tj =0, j = 2,3,...,s, (12) 



A,_i 



and 



n 9 



y-i 



A,_i 2 . 



-J>-i + l)4 = 0, J = a + l,...,fc- (13) 



Let us ?/y define as 

Vij = {k-j + 1) *«%«] + 2 ( fc - •?' + 1)4%^]' ( 14 ) 

where J^ is an indicator function whose value is one if A is satisfied and is zero if A 
is not satisfied. Then (ITTj) . ( TT2l and ffTB"]) can be expressed as 



, n k 

y-££Viv* = o, (is) 

i=i j=i 



and 



Vi = "T- ~ ^ £ !/y = °. j = 2, 3, . . . , fc. (16) 

Notice that, by using (114J1 . we have combined the two equations in ( 112 j) and (fT3l) 
which resulted in the equation (fT6|) . 



7 



For convenience, let y.j = Y^i=\Vij- Using this and considering Ao 
rewrite (IT5|) and ( JT6l) as follows: 



1, we can 



nk 

T 



y 



•1 ~ 5Z A i-!^' = ' 



i=2 



and 



(17) 



/? 



b-i 



A,_! 



%.j = 0, i = 2, 3, . . . , k. 



Solving f fT8|) for A,_i, we have 



Aj_i 



/? 



Oy.j 



j = 2,3,...,fc. 



(19) 



Substituting (!T9|) into (TT71) gives 



nk 

-j-y.i 



n(k- I) 




(20) 



Solving ( 1201) for 6, we obtain the MLE of 8, denoted by 6, 



n 



n 



En 
t=l 2/il 

The MLEs of Aj_i are also easily obtained by substituting (I2~T]) into (TT9|) 

En 
»=1 ^ x „• _ o q j, 

Aj_l — — _ n , J — Z,, O, . . . , fv. 

y»j i^i=i Vij 



(21) 



Notice that given the definition of yij, we can rewrite y,j in the following manner: 

n n 

J2(k ~ J + l)Uj = {k-j + 1) J^, J = 2, 3, . . . , s 



8=1 



8=1 



- J> - j + \)t% = -{k - j + 1) X)t£, j = 8 + 1, 



2 

*■ 8=1 



8=1 



Therefore we can write the closed form MLEs for 9 and A 7 _i as 



and 



n 



kj^t 



ii 



i=l 



fcj^*ii 



j=l 



(k-j + l)Y,t 



'J 



Vl = < 



j=l 



t=i 



^-i+i)E*? 



«=i 



j = 2, 3, . . . , s 



j = s + l,...,fc 
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